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Abstract. We present early results of our multiwavelength study of the star form- 
ing region IRAS 20050+2720. While we use X-rays and IR to classify young stellar 
objects, the optical data can be used to exclude foreground objects. The dataset set con- 
tains 57 class I sources, 183 class II sources and 183 X-ray sources, of which 140 are 
class III candidates. Within IRAS 20050+2720 four subclusters are found. Subcluster 
A and B are the densest regions, which contain stars of all evolutionary stages. Sub- 
cluster C is much younger than the other two. It has not formed any class III objects 
yet. We newly identify a fourth subcluster, which consists mostly of class II objects and 
is located about 10' from the center of the cloud. 



1. Introduction 

Star formation occurs pre-dominantly in dense and cool molecular clouds, which col- 
lapse under their own gravity and fragment into smaller pieces. In our galaxy we ob- 
serve those clouds over a range of masses, but it now seems that the rtiajority of stars 
form in clusters with m ore than a hundred members ( Carpente"^ llOOOl : [Lada & Lada 



20031 : iPorras et"al]|2003l ). Thus, studying star formation in those clusters is an impor- 



tant step to understand the history of the stars we currently observe. 

In the early stages the proto-star is still hidden by its parent molecular cloud and 
thus can only be observed as far-IR emission from the warm, dusty envelope (class 0). 
Due to the conservation of angular momentum the matter does not collapse radially 
onto the star, but forms an accretion disk in class I objects. They often drive powerful 
outflows and carve holes in their envelopes so radiation can escape. Eventually, the 
envelopes disperse and the stars become visible as classical T Tauri stars (CTTS), or 
class II sources in IR classification. Their disks still cause an IR excess, and on the 
stars accretion shocks and coronal activity lead to X-ray emission. Later, the IR excess 
vanishes, because the disk mass decreases. In this stage (class III or weak-line T Tauri 
stars = WTTS) cluster members cannot be distinguished from main-sequence (MS) 
stars by IR observations; one method to identify them is through their X-ray luminosity 



stars by IK observations; one mettioa to laentiiy ttiem is ttirougti ttieir A-ray luminosity , 
which is far higher than for older stars (for a review see lFeigelson & Montmerldl 19991) . 



I 



2 



H. M. Gunther, S. J. Wolk et al. 



In dense clusters, stars may influence one another during their evolution. Most 
notably strong winds or strong radiation fields of early-type stars can evaporate the 
disks of their late-type neighbors. In the Orio n nebula this proces s can be directly 
observed in distorted disks, called proplyds (e.g. lO'dell & Wen[|l994h . 

Despite significant observational progress many details of the star formation pro- 
cess in clusters are still under debate. The variety of competing processes makes it 
difficult to quantify the different contributions. For example, it is unclear if early or 
late-type stars form first in an undisturbed cluster. The theoretical assumption of a 
cluster evolving with r ninimal influenc e from the rest of the galaxy may be far from re- 
alistic. For e.g. p Oph. lde GeusI ( 1992 ) suggested that the expanding shell of the upper 
Sco star forming region traveled through the cloud and triggered star formation. If it 
turns out that this is the rule rather than the exception, it makes it even more difficult to 
disentangle the competing processes that influence star formation. 

In this situation, we decided to observe IRA S 20050-1-2720, wh ich is a star form- 
ing region located at a distance of about 700 pc (Wilkin g et al.lll989i) . Apparently, no 
massive star formed in IRAS 20050-1-2720, thus the intensity of the ambient radiation 
should be small and we can study the evolution of late-type stars in the absence of 
external irradiation on the disk. Observationally, this region is especially suited for 
our study, because the diffuse background at 24 /urn is weak, thus Spitzer observations 
should provide a more complete source list than in other clusters. 

IRAS 200504-2720 was initially discovered as a point source by the IRAS satellite. 
Molecular line e mission indicates a mass infall on the regions (Gregersen et al. 19971; 



Choi et al.lll999f). The lumi n osity of this region is estimated to 388 Lq from the IRAS 



fluxes ( Molinari et al.|[l996h . IChini et aTl ( 2001 ) obtained mm-maps and they estimate 
the total gas mass within 65"=0.2 pc radius as 200 Mq. 

The first attempt to establish an IR clas sification lead to about 100 cluster mem- 
bers, about half of which show an IR excess ("Chen et alJli997|). They find four regions 
of enhanced stellar density and label the most significant ones as A, B and C. Subcluster 



A con tains several deeply embedded YSOs (young stellar objects). While IChen et al. 



(119971) estimate an average cluster age of 1 Myr, they interpret the presence of radio 
lobes and the stronger reddening as signatures for a more recent star formation event in 
subcluster A in the last 0. 1 Myr. iGutermuth et all (l2005l) confirmed these ideas, noting 
that subcluster B lacks the 850 //m emission which is present in A and C. It seems, that 
subcluster B recently cleared the gas of the parental molecular cloud. 

With Spitzer we can get very accurate photo metry o f star forming regions and clas- 
sify individual objects according to their IR SED. lGutermuth et al.l(i2009() identified 177 
YSOs in IRAS 200504-2720, which could be grouped in two distinct cores. We extend 
the spatial coverage with additional Spitzer observations and add information from a 
deep Chandra observation and, surprisingly, the first dedicated optical photometry. 



2. Observations 



2.1. Spitzer 



Spitzer has observed IRAS 20050-1-2720 several times. We use the data o f two initial 



pointings (AOR ID: 3656448 and 3665 152), which are described in detail in lGutermuth et al 



(120091). Because this dataset covers only parts of the cluster follow-up observations 
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were performed to extend the coverage of Spitzer data with AOR 21893120 and 21893376 
(IRAC) and AOR 3665152 and 21891840 (MIPS). 

All datasets were reduced and combined with the same procedures as in j 
(12009) using the ClusterGrinder IDL package, which is described in detail in|( 



Gutermuth et al 



Gutermuth et al 



(12008). The procedure matches IRAC, 2MASS and MIPS sources and delivers a com- 
bined source catalog. The IR SEDs are classified using all available mid-IR color in- 
formation. 

2.2. Chandra 

IRAS 20050-1-2720 was observed with Chandra/ AClS-l. The total exposure time of 
95 ks was split in three exposures (ObsIDs: 6438, 7254 and 8492) with individual 
exposure times of 23 ks, 2 1 ks and 5 1 ks, respectively. The observations were taken on 
2006-12-10, 2006-01-07 and 2007-01-29. 

The individ ual exposures were me rged and processed with the pipeline of the AN- 



CHORS project (iSpitzbart et al.l 120051) . which aims to process all young star forming 
regions in the Chandra archive in a consistent manner Q. In short, sources are detected 
with the wavdetect algorithm in CIAO, and for each source event lists and response 
files are generated. For all sources with more than 20 counts a single Raymond-Smith 
thermal emission model is fit, which contain as free parameters the absorbing column 
density Nh, the temperature T and the volume emission meassure. The metalicity is 
fixed at 0.3 compared to solar. In addition, for bright sources with more than 100 counts 
a model with two temperature components is fit. Also, for ea ch source a lightcurve is 
generated and its Gregory-Loredo statistic is calculated dGregory & Loredo,,1992) to 
characterize the source variability. 

In total 1 83 sources are detected with Chandra. 

2.3. Optical photometry 

Optical photometry of IRAS 20050-1-2720 was performed on 2009-1 1-I8 from the Fred 
Lawrence Whipple Observatory (FLWO) with the 1.2 m telescope and the Keplercam 
instrument. Keplercam is a single-chip CCD camera with a Fairchild CCD 486 detector 
with a total field-of-view (FOV) 23.1 x 23.1'. Each pixel corresponds to 0.336"on the 
sky. We operated the chip in the default 2x2 binning. Observations were taken with 
staggered exposure times of 6, 60 and 300 s, where photometry for brighter stars can 
be performed on frames with shorter exposure times. Keplercam is equipped with a set 
of filters, here we use observations taken in BuVuri. The Bu and Vh filters belong to 
the Harris set, the r and / filters are part of an SDSS set. 

We performed the data reduction in Pyraf, a Python inter face to IRAF (lTodyll993h . 



Photometry is done with PSF-fitting in DAOPHOT (Stetson] [19871) as implemented in 



IRAF. Data and flatfield are bias-corrected, then the data images are divided by the 
median-averaged flatfield. The tail of the PSF is very wide. Thus, our source lists are 
incomplete close to bright sources and the completeness limit is non-uniform over the 
FOV. 

We correct the WCS of each exposure with the IRAF task maccmatch, which min- 
imizes the residuals between sources in our images and the 2MASS catalog. Sources 
are extracted independently for different exposures and source lists are merged with 
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a matching tolerance of 1". Sources, which are identified in exposures with different 
exposure times, are taken from the shorter exposure. 

We attempt to perform absolute photometry. Extinction and zero point o f the in- 



sttumental magnitude were calibrated with observations of Landold field 01 (iLandoh 
Il992h . which was reduced analogues to the data of IRAS 20050+2720 itself. We do not 
perform explicit conversions between the Harris (Bh and Vh fil ters) and the SDSS sys- 



tem (r and / filters) and the color system used by lLandoltl (119921) . instead, we just fit the 



instrumental magnitudes in BuVnri to the catalog values for BVRI assuming an offset, 
a linear color term and a linear extinction term. The color term absorbs the linear color 
transformation between the different systems. We tested this simplified procedure on 
stars in the standard field and we find that over the range of colors present in the Lan- 
dolt standard stars, all residual color trends are negligible compared to the uncertainties 
resulting from the non-uniform background and the PSF fitting. We remove the faintest 
stars from our final sourcelist, so that the typical error on the absolute flux is about 10%. 

2.4. Source matching 

The different infrared observations are matched as in Gutermuth et al. I (120081) . We then 



add the X-ray data, matching X-ray sources and IR sources, if they are closer than the 
large axis of the error ellipse from wavedetect. If more then one IR source fullfills this 
criterion, we match the closest X-ray and IR-source. Optical sources are identified with 
IR or X-ray sources using a constant matching radius of 3". The dataset contains 57 
class I sources, about 10% of them (6) are detected in X-rays, and 183 class II sources, 
20% (41) of them are detected in X-rays. The detection fraction of class II sources is 
higher, because more evolved sources have less circumstellar absorption. 



3. Results 

3.1. Class III sources 

The main goal of this study is to extend the treatment of YSOs to more evolved class III 
objects. While class I and class II objects can be identified by their IR excess, class III 
object have lost their envelope and their disk already. In the IR they are indistinguish- 
able from MS-stars. However, because of their low age, they are still fast rotators and 
consequently luminous X-ray emitters. Those X-ray sources, which do not coincide 
with either class I or class II sources, thus provides a list of class III candidates. 

3.2. Optical data 

Most of the YSOs in IRAS 20050-^2720 are not detected in our optical dataset. The 
limiting magnitude is about 19 in the r-band, which is insufficient for late-type objects 
in 700 pc, if they are reddened by a few magnitudes of dust in the proto-stellar cloud 
(We could detect a 0.5 Mq object behind Ay = 5 mag, but only 3.5 Mq and earlier 
behind Ay - 10 mag). The detection efficiency is much better in the IR, where the dust 
opacity is lower. There is very little overlap between the optical sources and IR class II 
sources or even the class III candidates. We can make use of this fact to identify X- 
ray emitting foreground objects to remove ffiem from the sample of potential class III 
sources. In the RIR - I color-m agnitude diagram (Fig. [B the extinction moves stars in 
parallel to the 1 Myr isochron e ( Siess et al]|200Cll) . We also show a reddening vector for 



Ay = 5 dSchlegel et al.ll 19981) . The values corresponds to the age of IRAS 20050-1-2720 
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Figure 1 . Color-magnitude diagram of all optically detected sources 



dChen et alJll997h and a typical value for the thickness of the cloud (iGutermuth et al. 
20081). The isochrone matches the class II sources reasonably well, but there are about 
a dozen X-ray sources to the left of the isochrone in the diagram, which cannot be 
explained by reddening; sources in this region must be foreground objects, which can 
be confirmed by overplotting the main sequence for distances of e.g. 300 pc. After we 
remove all X-ray sources with optical colors R - I < 0.5 for R < 16 and R - I < 0.75 
for R > 16 from our list of class III candidates 140 objects remain. 



3.3. Distribution of sources 
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Figure 2. Class I, II and III sources (symbols as in Fig [3]). The background im- 
age is Spitzer/IRAC 3.6 fj.m, the white contours show Ay from a 2MASS extinction 
map in steps of 2 mag (IGutermuth et alj|2005b. Red boxes show the position of the 
subclusters A, B and C as defined bv lChen et al.l (il997.) and the new subcluster D. 
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Fig. |2] shows a map of IRAS 20050+2720, where class I, II and III sources are 
overlaid on a Spit zer/IRAC 3.6 /urn image. We identify the subclusters A, B, C and D. 



Chen et alj (Il997f) named the three subclusters A, B and C, where A and B are nearly 



continuous. These clusters have the highest density of sources. Nearly all class I and 
most class II sources clump here. Subcluster C seems much younger than the other 
regions of the cloud, because not a single class III candidate is found there, indicating 
a very recent star formation event, which has not left enough time for stars for evolve 
to the class III stage. Conversely, the newly identified subcluster D is older than C 
as it does not harbor class I sources any longer, on the other hand very few sources 
have advanced to class III either. Also, the column density of the cloud is lower than 
in the other subclusters, possibly because stellar feedback contributes to its dispersal. 
The class III candidates are more spread out and are also found in the region between 
subclusters ABC and D. They might have been formed by earlier star formation events 
or they were ejected from the denser regions over time. 



3.4. Disk properties 
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Figure 3. Infrared color-color diagram with data from 2MASS. Sources in the 
bottom right region have optically thick disks. 



An infrared color-color diagram, such as Fig. [3j allows us to search for optically 
thick disks because the reddening vector is largely parallel to the isochrone. Reddening 
due to the proto-stellar envelope or the diffuse could gas alone displaces objects only 
along this vector, sources in the lower right region of the diagram must therefore have an 
excess in H and K-band emission compared to the other members of the sample. Extra 
emission at these wavelengths is caused by an optically thick inner disk. This region of 
the diagram is, as expected, populated mainly by class II sources because most class III 
sources have lost their disks already or at least their disks are vastly depleted. 



3.5. Reddening and extinction indicators 

In this section we attempt t o compare the thickness of the cloud, expressed as optical 
extinction Ay map ( Guterm uth et al.ll2005h . with the extinction of individual stars. For 
stars without circumstellar material we should find a distribution of individual Ay val- 
ues between and Ay of the cloud at that position, reflecting the distribution of sources 
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along the line-of-sight. However, the real situation is significantly more complicated. 
Sources of error include: Young stars, even in class III, still have circumstellar material, 
which gives them extra absorption; the spectral energy distribution has photospheric, 
accretion and disk contributions, which makes is difficult to measure the reddening; and 
the Ay map is made from averaging different sources including potential cluster mem- 
bers, which may be located in front of the cloud or have circumstellar material. All this 
introduces a large error on the values for individual stars shown in Fig. ID From a sta- 
tistical point of view it is reassuring, that most stars cluster around the blue line, which 
indicates equal values for both the reddening of the star and the thickness of the cloud. 
The scatter is large, but in general more absorbed objects are located in regions where 
the cloud is thicker, consistent with figure ^ The measurement is not good enough to 
trace the distribution along the line-of-sight. 



25 




Ay of source 



Figure 4. The reddening towards class I, II and III sources compared to the col- 
umn density of the cloud at that position. The blue line indicates equal values for 
cloud thickness and individual absorption. 



4. Summary 

We presented data on the young stellar cluster IRAS 20050-1-2720, which includes X- 
ray data from Chandra, optical data and IR data from 2MASS and Spitzer. The optical 
data is used to separate out foreground X-ray sources. Class I and class II sources are 
identified from their IR spectral energy distribution (SED). We treat all X-ray sources 
which are neither foreground objects nor classified as class I or class II according to 
their SED as potential class III objects. The spatial distribution shows a strong sub- 
clustering with four subclusters. Subcluster A and B have the highest density of young 
stellar objects with members of all three classes. Subcluster C is so young that no 
class III objects are present. Subcluster D is separated from the other three by about 
10'(=2 pc) and contains almost exclusively class II objects. We identified some objects 
with optically thick disks. 
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